Model of C-Axis Resistivity of High-T c Cuprates 



Yuyao Zha, S. L. Cooper and David Pines 

Department of Physics, University of Illinois at Urbana- Champaign, 1110 West Green Street, 

Urbana, IL 61801 
(February 1, 2008) 

Abstract 

We propose a simple model which accounts for the major features and 
systematics of experiments on the c-axis resistivity, p c , for La2_ x Sr x Cu04, 
YBa2Cu306+ x and I^S^CaC^Os- We argue that the c-axis resistivity can 
be separated into contributions from in-plane dephasing and the c-axis "bar- 
rier" scattering processes, with the low temperature semiconductor-like be- 
havior of p c arising from the suppression of the in-plane density of states 
measured by in-plane magnetic Knight shift experiments. We report on pre- 
dictions for p c in impurity-doped YBa2Cu3C>6+x materials. 
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Although there is currently no consensus ]TJ as to the important mechanisms contributing 
to c-axis transport in high temperature superconductors, recent transport and optical exper- 
iments of Bi2Sr2CaCu20 x , YBa2Cu306+ x , and La2_ x Sr x Cu04 |2|-|J reveal a number of key 
features that must be accounted for in any successful model of the interlayer charge dynam- 
ics in the layered cuprates: First, p c {T) in E^S^CaC^Og and underdoped La2_ x Sr x Cu0 4 
and YBa2Cus06+ x have a semiconductor-like temperature dependence (dp c /dT < 0) at low 
temperatures, and a linear-in-T dependence at high temperatures. The crossover tempera- 
ture, T*, between these two regimes decreases with increased doping in YBa2Cus06+ x and 
La2_ x Sr x Cu04. Importantly, c-axis optical measurements || show that the semiconductor- 
like resistivity "upturn" in underdoped YBa 2 Cu30 6+x is actually associated with a uniform 
suppression of the optical conductivity below ~ 300cm -1 . These data suggest that the c- 
axis conductivity scales at low frequency with the Knight shift, which is proportional to the 
in-plane density of states: K s oc N(0). Second, both La2_ x Sr x Cu04 and YBa2Cus06+ x ex- 
hibit a strongly doping-dependent mass anisotropy, and a possible crossover from quasi-2D 
to 3D transport behavior at high doping, that arises in part from doping-induced structural 
changes [§,[7J- 

In this paper, we present a simple phenomenological model of the c-axis resistivity in 
the layered cuprates that describes the key elements of interlayer transport in the cuprates 
with reasonable parameters. Before examining the mechanisms contributing to c-axis trans- 
port in the cuprates, it is important to point out that experiment evidence suggests that, 
except perhaps for the overdoped cuprates, c-axis transport in high T c superconductors is 
incoherent. For example, typical estimates of the c-axis scattering rate in YBa2Cu306+ x 
give 1/t c > 1000cm -1 0, while optical measurements indicate a c-axis plasma frequency 
of u p ± ~ 40meV in fully oxygenated YBa 2 Cu 3 07, and u p ± < lOmeV in underdoped 
YBa 2 Cu30 6+x |7J]. These values suggest that the c-axis mean free path is of order or less than 
the c-axis lattice spacing, i.e., c-axis conductivity in the cuprates is below the Ioffe-Regel 
limit, and hence c-axis transport is incoherent. 

One important contribution to c-axis transport in the cuprates is expected to arise from 
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electron scattering in the "barrier" layer between Cu02 "cells" (i.e., layers, bilayers, trilayers, 
etc.). For example, c-axis Raman scattering measurements provide evidence that carriers 
hopping between layers scatter from c-axis optical phonons associated with the barrier in 
YBa2Cu306+ x 0j while Littlewood and Varma have pointed out the likelihood that static 
impurities in the barrier layer provide an important source of scattering for c-axis transport 
in the cuprates [|TJ . A phenomenological expression for this contribution to c-axis transport 
can be written as: 

a® = N(0)^ftlr c (1) 

where d is the interlayer spacing, N(0) is the in-plane density of states, t± is the interlayer 
coupling, and r c is the c-axis scattering time. 

On the other hand, c-axis transport measurements of La2_ x Sr x Cu04 and Bi 2 Sr 2 CuO x 
yield p c oc p a b at high temperatures [||,|J, suggesting that scattering or fluctuations in the 
planes may dominate c-axis transport in this regime. This contribution to c-axis transport 
can be written [[J, 

= N(0) e ^tjr ab (2) 
where r a b can be derived from the planar conductivity, 
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and the temperature-independent quantity, t±, measures the effectiveness of planar scattering 
processes to c-axis transport. 

Because Eqs.([L]) and (Q) describe independent physical processes, it seems natural to 
consider the corresponding tunneling and/or scattering mechanisms as additive in the resis- 
tivity. We are thus led to propose the following expression for p c , 

h 2 1 1 

Pc N(0)e^ [ tlr ab + tlrj- W 

In the limit that one mechanism or the other is dominant, Eq.(^) yields the corresponding 
conductivity given in Eqs.(|I]) and (Q). We now proceed to use Eq.([|) to analyze c-axis 
transport experiments. 



We consider first the La2_ x Sr x Cu04 system. We obtain the temperature-dependent den- 
sity of states N(0), from the recent analysis by Barzykin, Pines and Thelen [PH|1 , who extract 
the temperature-dependent uniform susceptibility, Xo{T), from Knight shift measurements 
and scaling arguments on this system; in deriving N(0), we neglect Fermi liquid corrections, 
determining N(0) from Xo{T) = —fi%N(0) where ps is the Bohr magneton. The density 



of states N(0) obtained from the Knight shift data of Ohsugi et al JTTj] is shown in Fig. 
1. We determine r a b by using Eq.([|) and optical measurements of ui p \\ 0. An independent 
measure of t±, the interlayer hopping amplitude, can be obtained from the c-axis plasma 
frequency, u> p ±, measured in optical or penetration depth experiments J7|. This leaves us 
with a two-parameter fit to the data for each hole doping, these parameters being t\ and 
l/r c . Our results are given in Fig. 2 and Table I, and we comment on them briefly. 

We first note that at each doping level, the high temperature behavior of p c is deter- 
mined entirely by the planar conductivity, according to Eq.(@), so that t\ is completely 
determined by p^. The values of t\ obtained in this way are independent of hole concen- 
tration, within 10%. Secondly, both t± and r c are independent of temperature, and the 
barrier layer scattering contribution described by Eq.(|I|) dominates at sufficiently low tem- 
peratures. As might have been expected, r c is essentially independent of doping, while t±, 
which increases with increasing hole concentration, displays a strong dependence on hole 
doping where t± ~ x a with a > 2. According to Eq.(|4]), the crossover temperature T*, 
which separates the dp c / dT < and dp c / dT > regimes, occurs when the two terms in the 
parenthesis of Eq.(^) are equal. After some simple algebra, we find that at T*, we have 

p ab (T*) N T *{Q)eWHlTjl K u>l f J 2vriV T «(0)e 2 ^ {> 

where A^*(0) is the density of states at T=T*. Roughly, both A^»(0) and u p \\ increase 
with doping concentration, while all the other quantities on the right hand side of Eq.([5]) 
are doping independent. It turns out that the anisotropy ratio p c j p a b at the crossover 
temperature T* is nearly independent of hole doping for each individual cuprate. Thus 
for La 2 _xSr x Cu0 4 , Nakamura et al || found p c / Pat ~ 300 at T=T*, for doping levels of 



x = 0.10, x = 0.12, x = 0.15, and x = 0.20. 

We now apply the same method of data analysis to YBa2Cu306+ x experiments. Our 
results are shown in Figs. 3 and 4, and Table I. In Fig. 3, we give our results for N(0), 



obtained from the analysis by Ref. [H| of the Knight shift experiments of Ref. ||12|| . The 
nearly doping independent results for r c shown in Table I are taken from fits to optical and far 
infra-red experiments || . Again we see that t\ is nearly independent of hole concentration, 
while t± increases with hole doping as t± ~ x a where a > 2. The anisotropy ratio at the 
crossover temperature T* remains almost doping independent. According to the data of Ref. 
@, this ratio is p c j p a b ~ 100 at T=T* for different doping levels. 

A further test of our model is provided by the c-axis resistivity measurements on 
Bi2Sr 2 CaCu 2 08 shown in Fig. 5, which display the same qualitative dependence on tem- 
perature. Our theoretical results are obtained by deducing N(0) from the Knight shift 



experiments of Walstedt et al |14| , r a b from the measurements by Martin et al of p a b and t± 
from Ref. 0. The agreement between our model calculation and experiment is again seen 
to be satisfactory. The parameters used in making the fit given in Table I. 

It is instructive to compare the variations in ti and t±_ on going from one system to 
another. The fact that t\ is the largest in YBa 2 Cu306+ x , somewhat smaller in La2_ x Sr x Cu04 
and smallest in Bi 2 Sr 2 CaCu208 suggests that although independent of doping, the dephasing 
of the afe-plane enters p c in a way which is related to the interlayer (unit cell) distance in 
these cuprates. Another observation from Table I is the systematic behavior of the c-axis 
"barrier" scattering rate l/r c : the scattering is strongest in YBa 2 Cu 3 6+x , which contains 
two apical oxygens per unit cell; weaker in La 2 _ x Sr x Cu04, which has one apical oxygen per 
unit cell; and significantly reduced in Bi 2 Sr 2 CaCu 2 08, which does not contain any apical 
oxygens. Our result for l/r c is also consistent with the Drude fitted results from optical and 
Raman experiments for YBa 2 Cu 3 6+x All these facts support the suggestion that the 
c-axis barrier scattering originates most likely from the c-axis apical oxygen phonons. 

c-axis transport measurements on impurity-doped systems provide a direct test of our 
description of p c in Eq.(f|). We consider two types of impurities: Zn and Co, both of which 



have been doped into YBa 2 Cu307. Zn is known to go in as a planar substitute for Cu. To first 
order, it does not modify the planar density of states N(0), or the in-plane scattering rate 
l/r a b, nor will it influence t± and r c . Therefore, we predict that for Zn doped YBa 2 Cu307, 
the changes which take place in the c-axis transport will mirror the comparatively small 
increase, 5p a b oc nz n , found for this system JTH]. On the other hand, when Co is doped into 



YBa 2 Cu307 up to a 2.5% doping level, it does not influence T c [pT6| , and may plausibly be 
assumed to substitute for chain Cu atoms. It will therefore not affect t±, T a b, or N(0), and 
we may expect Eq.(@) to be unaffected by Co substitutes to this level. Nevertheless, 
we would expect that t±, which is sensitive to the chains, will be reduced, thus increases 
the magnitude of the barrier scattering contribution to p c , and we expect p c to have a 
temperature dependence somewhat similar to that of YBa2Cu307: linear-in-T for the whole 
temperature range, with the same slope but a larger residual resistivity at T=0. 

In conclusion, we have presented a model for c-axis resistivity of high T c cuprates, in 
which both in-plane dephasing and the c-axis "barrier" scattering contribute to the c-axis 
resistivity. The behavior of the "barrier" scattering is consistent with that of the c-axis 
apical oxygen phonons. Our model fits quite well with the existing data on La 2 _ x Sr x Cu0 4 , 
YBa 2 Cu 3 6+x and Bi 2 Sr 2 CaCu 2 8 . 

We acknowledge stimulating conversations with A. J. Leggett, Q. M. Si, A. G. Rojo, D. 
M . Ginsberg, K. Levin, J. T. Kim and L. Y. Shieh. This work is supported by NSF-DMR- 
91-20000. 
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FIGURES 

FIG. 1. The planar density of states of La2- x Sr x Cu04 obtained from the Knight shifts data 



of Ohsugi et al [11|, following the scaling analysis of Ref. 



FIG. 2. Calculated p c for La2- x Sr x Cu04 at different doping levels (symbols), plotted against 
the experimental data of Nakamura et al Q (solid lines). 

FIG. 3. The planar density of states obtained from the analysis of Ref. |l3| for YBa2Cus06+ x , 
here we use N(0) of Yo.gPro.iBaaCusOr for YBa 2 Cu 3 6 . 88 . The N(0) for YBa 2 Cu 3 6 . 68 is an 
estimate from the scaling arguments given in Ref. fl3|| . 

FIG. 4. Calculated p c for YBa2Cus06+ x at different doping levels (symbols), plotted against 
the experimental data of Takenaka et al Q (solid lines). 

FIG. 5. Calculated p c for Bi2Sr2CaCu20s (circles) in comparison with experimental data by 
different groups 0. iV(0) is obtained from the scaling analysis of the Knight shift measurement of 
Ref. g4|. 
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TABLES 
TABLE I. Relevant Numbers 





Wp|| 


h 


V r a6 


t± 


1/Tc 


(t\T ah y l [at 300K] 






[eV] 


[meV] 


[meV] 


[meV] 


[meV] 


[meV- 1 ] 


[meV- 1 ] 


YBa 2 Cu0 6 . 68 


0.8 


14.4 


~ 2kT 


7.3 


100 


0.25 


1.88 


YBa 2 Cu0 6 . 78 


1.0 


14.1 


~ 2A;T 


10.4 


100 


0.26 


0.92 


YBa 2 Cu0 6 .88 


1.18 


18 


~ 2kT 


16.7 


100 


0.16 


0.36 


YBa 2 Cu0 6 .93 


1.4 


16.6 


~ 2£;T 


30-40 


109 


0.19 


0.089 


Lai.9oSr .ioCu0 4 


0.44 


3.37 


~ kT 


0.9 


6.78 


2.29 


8.37 


Lai. 88 Sr .i2CuO 4 


0.57 


2.78 


~ kT 


1.3 


6.99 


3.36 


4.14 


Lai. 85 Sro.i5Cu0 4 


0.7 


3.11 


~ fcT 


2.4 


6.56 


2.69 


1.14 


Lai. 80 Sro. 2 oCu0 4 


0.87 


3.01 


~ fcT 


5.5 


6.67 


2.87 


0.22 


Bi 2 Sr 2 CaCu 2 8 


1.25 


0.88 


~ kT/2 


0.1 


2.64 


16.8 


264 
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